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Abstract 

By introducing an additional operator into the action and using the 
Feynman-Hellmann theorem we describe a method to determine both the 
quark line connected and disconnected terms of matrix elements. As an 
illustration of the method we calculate the gluon contribution (chromo- 
electric and chromo-magnetic components) to the nucleon mass. 
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1 Introduction 



One of the earliest experimental indications that the nucleon consists not only of 
three quarks, but also has a gluonic contribution came from the measurement of 
the fraction of the nucleon momentum carried by the quarks. That this did not 
sum up to 1 as is required from the energy-momentum sum rule gave evidence 
for the existence of the gluon. Denoting (x)f as the fraction of the nucleon 
momentum carried by parton / we have 

X>>« + = 1 ' a) 

q 

where for the quarks / = q — u,d, . . . and for the gluon / = g. Experimentally 
(x)u+d ~ 0.4 so the missing component is large ~ 50% of the total nucleon 
momentum. Both (x) q and (x) g have similar definitions and so analogously to 
the definition of (x) q we have, with M. denoting Minkowski space 

(N(p)\[d M ^ lfi2 - \r)^6 M{9)a a \\N{p)) = 2(x) g [p^p^ 2 - \^m 2 N ] , (2) 
where 

(where 0{t) = f d 3 xO(t,x) and with normalisation (N(p)\N(j>')) = 2E N S(p — 
p'))- Note that we can generalise from a nucleon to an arbitrary hadron (averaging 
over polarisations if necessary). Higher moments can also be considered, by 
inserting covariant derivatives between the Fs. These occur when using the 
Wilson operator product expansion which relates them to moments of structure 
functions in a twist expansion. 

There have been many lattice estimates of the quark momentum fraction (x) q 
both for the nucleon (see e.g. [U [2] for a review) and the pion e.g. [31 H], but few 
attempts for the gluon part, (x) g EJ [7]. This is due to the fact that a lattice 
simulation must compute a quark line disconnected term, which is extremely noisy 
and gives a poor signal. These are direct calculations; in this letter we propose a 
new method using the Feynman-Hellmann theorem, to determine the gradient of 
En as a function of a parameter of an operator which has been introduced into 
the action S — > S(X) = S + XSo- An obvious disadvantage of this method is that 
it requires dedicated simulations for each operator of interest, but the gain, as 
we shall see, is a much cleaner signal. 

While the method is general, we shall demonstrate its practicability here by 
determining (x) g in the quenched case. 

2 The Feynman-Hellmann theorem 

We first briefly describe the Feynman-Hellmann theorem, in a Euclidean form 
that will be useful for the case to be considered here. Let S depend on some 
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parameter A, so S — » 5(A). Now as by definition the (Euclidean) correlation 
function is given by 



(the unpolarised case for the nucleon and where we make the obvious replacements 
N by H and N by for other hadrons), then we have 

We now use the transfer matrix formalism on both sides of this equation. Ignoring 
finite size effects this gives 

(N(t)N(0)) x = A N (X)e- EN{x)t + exp. smaller terms . (6) 

so on the LHS of eq. 

■^-(N{t)N{0))x = - (N(t)N(0))x t + exp. smaller terms. (7) 

OA OA 

Furthermore, if f2(r) is any operator (local in time), then using the transfer 
matrix formalism again the associated 3-point function gives 



(N(t)n(r)N(0)) x _ { 2E^X)( N P\ N )* + ex P- sma11 terms < T < 1 
(N(t)N(0)) a I exp. small terms otherwise 



(8) 



Note that we have inserted a 2En in the denominator of the RHS to account 
for the mis- match of normalisations, i.e. to agree with those of eq. ([2]). Hence 
summing over r also gives a linear term in t. Thus from this equation, replacing 
^(t) by the operator in the RHS of eq. (jSJ), and together with eq. ([7j) we have 
the Feynman-Hellmann theorem 



dE N (X) 1 , 

d\ 2E N {X) 



dSjX) 
dX 



N , (9) 



(where : . . . : means that the vacuum term has been subtracted). Thus by suitably 
choosing So and by identifying numerically the gradient of E N (X) at A = we 
can determine the desired matrix element. 
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3 The lattice method 



3.1 Gluon operators 

Before considering the lattice, let us first Euclideanise the gluon operators^ to 
give us an indication of what we might add to the action. Defining 

Ofiu = tr cF^iaFua i (10) 

(tr c F 2 = \F a2 ) this then gives the two obvious operator choices (a) and (b), 
O ai = O i4 = tr c (E x B)i 

O b = O u - \0 3j = §tr c (-£ 2 + B 2 ) (11) 

^qM(9) _^ ^ M{g) _^ ^ re l a tion to (x) g is given by 

(N(p)\d at \N(p)} = -2iE NPi (x) g 

(N(p)\d b \N(p)) = 2(m 2 N + lp 2 )(x) g , (12) 

with 

-> -> ->2 -»2 

O ai = tr c (£ x B)i , 6 b = §tr c (-£ +B). (13) 

Both choices have their difficulties: operator (a) always needs a non-zero mo- 
mentum p, while operator (b) requires a delicate subtraction between two terms 
similar in magnitude. 

Note that, because of Euclideanisation (footnoted]) the energy has a negative 
£ 2 term, while the action (see section I3.2p has a positive S 2 term. 

3.2 The action 

We now turn to the lattice. We shall use the Wilson gluonic action 

5=I/3^Retr c [l-f/°(a;)] , (14) 

(i.e. sum over plaquettes), with (3 = 6/g 2 . As 

Retr c [1 - U*(x)] = \a 4 g 2 F^(x) 2 + . . . , (15) 
J Our conventions follow [3]. So E Ml = F Mi0 -> iF u = iE { and B Ml = -\e ijk F M jk -> 
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this motivates the simplest definition of electric and magnetic field on each time 
slice as 

h£ a \r) = I/3l^Retr c [l-C/°(f,r)] 

CI 

xi 

\B a2 {r) = Rfitr e [l-0g(£,r)] , (16) 

x i<j 

respectively. For the action we thus take 

5(A) = a W a \r) + B a \r)] - Xa £ \[-£ a \r) + B a \r)] , (17) 

T T 

or in terms of the gauge plaquettes 

S(X) = ±/3(l + A)^Retr c [l-f/°(f,r)] 

i 

+i/3(l-A)^Retr c [l-C/g(x,r)] . (18) 

i<j 

Of course for A = 0, then this reduces to the standard action, eq. (fH|) . 
3.3 Gluon moment 

Comparing the results of sections 13.11 and 13.21 we see that they can be applied 
to operator (b) only; operator (a) would require the clover definition of the field 
strength tensor. Using eq. ([TT]) together with eq. (fT2|) and eq. ([9]) gives from the 
Feynman-Hellmann theorem 

which leads to 

dE N (X) 



dX 



A=0 2E N 



(m% + If) (x)? , (20) 



where the lat superscript on (x)™ signifies that it is now the lattice operator. 
The vacuum term which appears in section [2] has been dropped, because 

(0|±(-£ a2 + £ a2 )|0> = 0. (21) 

This follows from rotation symmetry. In the Euclidean vacuum the time and space 
directions are equivalent, so the average trace of the chro mo-electric plaquettes, 
U^, is the same as that of the chromo- magnetic plaquettes, C/R, in eq. ffTEj) . 
leading to perfect cancellation in eq. (l2Tj) . 
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4 Lattice results 



We work with quenched Wilson clover fermions at ft = 6.0, c sw = 1.769 and k = 
0.1320, 0.1324, 0.1333, 0.1338, 0.1342 on a 24 3 x 48 lattice with antiperiodic time 
boundary conditions for the fermion. We have generated O(500) configurations 
for each ensemble. We use standard nucleon interpolating operators together with 
Jacobi smeared source/sink as in e.g. [3]. The results were generated using the 
Chroma program suite, [SJ. We have only considered the case p = so eq. (1211 
reduces to 



dam at (A) 



3am a? 



(22) 



A=0 



To estimate the gradient at A = 0, we have generated data at A = —0.03333, 0.0, 



0.03333 which enables us to straddle the A 
given in Table [TJ 



point. The raw data results are 



K 


A = -0.03333 


A = 


A = 0.03333 


0.1320 


1.0033(29) 


0.9772(33) 


0.9564(34) 


0.1324 


0.9537(30) 


0.9283(34) 


0.9077(36) 


0.1333 


0.8357(33) 


0.8117(40) 


0.7923(41) 


0.1338 


0.7649(38) 


0.7413(47) 


0.7236(47) 


0.1342 


0.7044(47) 


0.6799(62) 


0.6647(55) 



Table 1: Nucleon masses, am^, as a function of A for five quark masses, k, calculated 
on ensembles with fixed ft = 6.0 and c sw = 1.769. 

In Fig. [l]we plot the nucleon mass, aniN, against A for the five quark masses. 
The data show no 0(A 2 ) effects for the A values chosen. These gradients (at 
A = 0) together with the nucleon masses (again at A = 0) determine (x) 1 ^ from 
eq. ( I22p which are given in Table [2j 



K 




{ X )T 


0.1320 


0.55499(48) 


0.4826(456) 


0.1324 


0.51745(49) 


0.4985(502) 


0.1333 


0.42531(52) 


0.5383(644) 


0.1338 


0.36711(55) 


0.5620(811) 


0.1342 


0.31433(62) 


0.5893(1062) 



Table 2: The pion mass and for the five different quark masses. 
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Figure 1: The nucleon mass against A for the five k values, together with a linear fit 
for each k value. 



5 Renormalisation 

As gluon operators are singlets, they can mix with the quark singlet. However 
there exists a combination of singlet operators with vanishing anomalous dimen- 
sion. (This is due to the conservation of the energy-momentum tensor, eq. ((!]).) 
We follow [B] and first write 



(x) b ; re + J2( x )T e = 1 + °( fl2 ) 



(23) 



where 



(x) b - = Z g (x) 



I a t 

9 ' 



{x)\- = Z q {x) l f . 



Together with the change to a scheme (here taken as MS) 



1 ~~ ZbareggW ^blreqqil 1 ) 



(x) bart 



(24) 



(25) 



this completes the renormalisation procedure. As we are considering quenched 
QCD only there is a simplification as Z^'J e = 1, 



(26) 
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{Z™J eqq {fi) is common for all the quarks). We thus need to determine Z g , Z q 
and Z^ eqq (fj,). We can find Z g by following [10] in considering an alternative 
interpretation of the action ( fl8j) . We motivated this action by adding a multiple 
of the gluon x operator to the standard action, but we could also write the action 

as 

S=lfrJ2 R etr c [1 - Uu&t)] + Retr c [l - U^(x,t)] . (27) 

i i<j 

which is the standard way of writing a gluon action on an anisotropic asymmetric 
lattice, with differing spatial and temporal lattice spacings, a s ^ at. This action 
has been studied in detail, in particular the way in which the anisotropy £ = 
a s /a t depends on (3 S and j3 t is known both perturbatively and non-perturbatively 
[TT] . At tree- level the anisotropy is given by £f ree = (3 t /(3 s - Z g can be found by 
comparing the anisotropy actually produced by splitting (3 S and f3 t with this tree- 

2 

level value. The result is Z g = 1 — y(c CT — c T ) where the anisotropy coefficients c a 
and c T are defined in [TT]. Using the perturbative values for c CTT [12] yields Z g = 
1 — 0.16677g 2 + • • • as the 1-loop perturbative Z g . In [9] this result was combined 
with non-perturbative determinations of c CTjT , [TT J , to give a Pade expression 

_ 1 - 1.0225^ + 0.1305^ 
Z °~ 1 - 07855V ' P ~ ' ( } 

(with an error of ~ 1%). So for /3 = 6.0 this gives Z 9 = 0.748. 

To estimate Z q we use the results for (x) l g at from Table [2] together with those 
for (x)^* from [T3] (i.e. u 2 b) together with eqs. (|231 and (|2TI) . In Fig. [2] 

we ploto (x)^ at + (x)^* against (x) l g at . From eq. ff23l we would expect that the 
y-intercept is given by 1/Z g and the x- intercept is given by l/Z g . At present 
we do not have enough results for a determination, so we shall just check for 
consistency by fixing the y-intercept as 1/0.748 and the x-intercept as 1, [B]. 
This gives consistency so we shall adopt here Z q — 1 together with a 10% error. 

Also from [13] , we have for [i = 2 GeV, 

Z^ reqq ^ = 2GeV)Z q = Z£x [AZ^(/i = 2 GeV)]' 1 

= 1.45 x 0.732(9) = 1.06(1) , (29) 

where the second equation uses the notation of that article (the non-perturbative 
RI — MOM scheme is converted to an RGI form and then back to the MS 
scheme). Further values of AZ^ S (^) are also given in [13]. With Z q this then 
gives Z^ eqq . 



2 The total contribution to (x) q from sea quarks has the form Nf x (disconnected term). So, 
even though the disconnected loop term is itself non-zero, we do not need to consider it because 
its coefficient vanishes if we work consistently in the quenched approximation. 
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Figure 2: 

y = (l - x)/0.748 



against (a;)" 1 * for the five k values, together with the line 



6 Results and conclusion 



We are now in a position to determine (x)" s (/i = 2 GeV). Using the first equation 
in eq. f l26|) together with eq. ( 128]) (evaluated at (3 = 6.0) and eq. (1291 gives 
(x)fS([x = 2 GeV). In Fig. [| we plot using eq. fl26]), (x)f^(fi = 2 GeV) versus 
(am^) 2 . This gives a value for (x)^ IS (fj, = 2 GeV) of 



(x)* /s (/i = 2GeV) = 0.43 (7) (5) 



(30) 



as our final result, where the first error is in the determination of (x) 1 ^ and the 
second is due to the renomalisation procedure. This is a significant improvement 
of our previous estimate 0.53(23) based on generating 0(5000) configurations, [5] 
(with error given just for (x) 1 ^). 

Direct measurements of gluonic expectation values are notoriously plagued 
by noise problems, because the gluons are bosonic fields. We have seen here 
that a cheaper alternative, modifying the gluon action and using the Feynman- 
Hellmann theorem to find expectation values from mass measurements, works 
well. Here we have performed a test calculation in the quenched case. The 
method is a generalisation of that used to determine the sigma term (see e.g. 
[T4"] and references therein), /3-function, e.g. |15] . or singlet terms, e.g. |l6j . It is 
clearly interesting to repeat this with dynamical fermions. 
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' lit H 
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Figure 3: (x)g IS (fi = 2 GeV) versus (om ff ) 2 for the five n values, together with a linear 
chiral extrapolation. 
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